Magnetic Resonance in Medicine

Purpose: Quantitative susceptibility mapping (QSM) provides venous contrast through the paramagnetic susceptibility of deoxyhemoglobin and can be used to estimate oxygen extraction fraction (OEF), a marker of cerebral metabolism. However, cerebral vein segmentation remains challenging due to artifacts, variability across QSM reconstruction methods, and limited sensitivity of conventional vessel-filtering approaches to small cortical veins. This study proposes a physics- and morphology-constrained deep learning framework for cerebral vein segmentation on QSM. Methods: Thirty subjects were manually segmented and used to train an attention-gated UNet. In addition to supervised Dice and cross-entropy losses, the network incorporated two selfsupervised constraints. A physics-informed loss enforced consistency between the measured local field and a field simulated from predicted veins using dipole convolution. A Frangi vesselness loss encouraged anatomically plausible tubular structures. Segmentation performance was evaluated on an independent multi-center cohort and compared with a prior vessel-filtering method and a supervised U-Net baseline. Physiological relevance was assessed by comparing QSM-derived venous OEF with calibrated fMRI-derived OEF. Results: The proposed model achieved superior segmentation performance (centreline Dice = 0.70 {+/-} 0.08) compared with MSVF (0.44 {+/-} 0.09) and the U-Net baseline (0.57 {+/-} 0.14). Improvements were most pronounced in smaller veins, resulting in greater vessel continuity. Venous OEF derived from the proposed framework demonstrated stronger agreement with calibrated fMRI-derived OEF than the previous approach. Conclusion: Combining supervised learning with physics- and morphology-based constraints improves vein segmentation and enhances the sensitivity of QSM-derived OEF measurements. Code and model weights are publicly available at https://github.com/QuantitativePhysiologyImagingLab/VeinSeg

Abstract: Purpose: Superparamagnetic iron oxide nanoparticles (SPION) are widely used as MRI contrast agents for cell tracking and molecular imaging. However, their detection using conventional T2/T2* weighted and susceptibility-weighted MRI relies on negative contrast (signal voids), which lacks specificity because similar hypointense signals can arise from other sources of magnetic susceptibility. This study introduces Frequency Shift Imaging (FSI), a positive-contrast technique based on an adiabatic Shinnar-Le Roux (SLR) self-refocused pulse that addresses the limited specificity of conventional SPION imaging. Methods: An adiabatic self-refocused RF pulse was designed using the SLR algorithm for frequency-selective positive-contrast imaging. The optimized pulse had a spectral bandwidth of 370 Hz, a peak RF amplitude of 14 T, and a duration of 15 ms. FSI data were acquired on a 7T MRI scanner. Phantom experiments used agarose-embedded 30 nm cell-labeled SPION at 0%, 25%, 75%, and 100% concentrations (pure). In vivo validation was performed following intramuscular SPION injection into the hind limb of a mouse. FSI images were compared with conventional acquisitions. Results: Quantitative phantom analysis showed normalized FSI signal intensities of 25.02% and 75.16% for samples containing 25% and 75% labeled cells, respectively, demonstrating an approximately proportional relationship between FSI signal and labeled-cell concentration. In vivo, FSI generated localized hyperintense positive contrast at the SPIO injection site, providing improved conspicuity and localization compared with the corresponding negative-contrast signal void. Discussion and Conclusion: FSI enables robust, B1-insensitive, positive-contrast SPION imaging at 7T with reduced echo time and SAR compared to paired adiabatic approaches, supporting its potential for quantitative cell tracking and molecular imaging applications.

Purpose: We aim to acquire and reconstruct high signal-to-noise ratio (SNR), whole-brain, mesoscale (isotropic submillimeter) diffusion MRI (dMRI) in vivo in a clinically feasible scan time. Methods: We present ROVER-dMRI, a unified rotating-view acquisition and super-resolution reconstruction framework. Thick slices (up to 10 mm)acquired across many rotated orientations improve the base SNR, and a shared multi-resolution hash-encoded implicit neural representation (INR) with per-direction MLP head and a physics-informed, Bloch-simulated slice-profile forward model reconstructs these views into a continuous isotropic volume (up to 13x super-resolution) in a self-supervised manner. ROVER-dMRI was evaluated in vivo on human brain (0.78 mm isotropic and 0.5 mm isotropic) and monkey brain (0.6 mm isotropic) across MRI platforms with varying gradient performance, and benchmarked against gSlider-SMS and least-squares super-resolution reconstruction (LS-SRR). Results: At 0.78 mm isotropic resolution, ROVER-dMRI diffusion- weighted images achieved higher SNR than the vendor product sequence, whose images were limited by the noise floor, while its directionally encoded color (DEC) maps showed lower noise than gSlider. When reducing the number of views from 20 to 13 (reducing the scan time by 65% ), ROVER-dMRI acquisition/reconstruction preserved all fine details in the fiber orientations. In the monkey brain, its DEC map matched a 60-direction gSlider acquisition at only 40% of the scan time. ROVER-dMRI preserved coherent fiber-orientation-distribution (FOD) structure, including crossing fibers, and recovered fine details in the cerebellum at 0.5 mm isotropic, that are not visible at 0.78 mm resolution. Compared to state-of- the-art 0.5 mm methods, ROVER-dMRI reduced scan time by 73% while acquiring more diffusion directions. Conclusion: ROVER-dMRI offers a practical route to high-SNR whole-brain mesoscale diffusion imaging in vivo, with the pulse sequence and reconstruction pipeline released as open source. KEYWORDS: diffusion MRI, super-resolution, implicit neural representation, hash encoding, physics informed learning

Abstract Purpose: To determine whether structural MRI information can improve reconstruction of highly accelerated 3D snapshot chemical exchange saturation transfer (CEST) MRI. Methods: Fully sampled brain CEST data were retrospectively undersampled at acceleration factors (AFs) of 4, 6, 8, 10, and 12. We compared reconstruction without structural information with population-level pretraining using T1-weighted (T1w), T2-weighted (T2w), or combined T1w-T2w images, as well as direct conditioning using a co-registered subject-specific T1w image. Performance was evaluated in two held-out healthy subjects using peak signal-to-noise ratio (PSNR), structural similarity index (SSIM), and mean absolute error (MAE) across frequency offsets. Z-spectrum fidelity and derived magnetization transfer ratio asymmetry (MTRasym) maps were additionally evaluated within three-dimensional gray- and white-matter masks. Results: Subject-specific T1w conditioning consistently achieved the highest PSNR and SSIM and the lowest MAE across AFs, with greater improvements at higher acceleration. It also improved Z-spectrum agreement, particularly in white matter and combined gray- and white-matter regions across most AFs. Improvements in MTRasym were more modest and varied across tissue types and acceleration factors. Conclusion: Incorporating subject-specific structural information can improve reconstruction of highly accelerated snapshot CEST MRI, particularly at higher acceleration factors. However, improvements in reconstructed source images did not consistently translate into better preservation of derived CEST contrast, highlighting the need for reconstruction strategies that more directly incorporate CEST spectral information.

The shift towards open science in preclinical research requires high-quality, comparable data adhering to FAIR principles, yet rigorous quality assurance (QA) and quality control (QC) frameworks remain less established in preclinical magnetic resonance imaging (MRI) than in clinical imaging. To address this limitation, a multicenter study was conducted across 21 international laboratories using standardized commercial liquid phantoms for mouse and rat MRI setups and a harmonized acquisition protocol. Data processing was centralized using AIDAqc, an automated pipeline extracting quantitative metrics including signal-to-noise ratio (SNR), temporal SNR (tSNR), ghost-to-signal ratio (GSR), and motion-equivalent temporal instability, combined with five complementary outlier-detection algorithms. The evaluated datasets encompassed magnetic field strengths from 3.0 to 16.4 T and heterogeneous coil and acquisition configurations. In mouse phantom data, SNR and tSNR were significantly lower at 3 T than at 7 and 9.4 T, whereas artifact-related metrics did not differ significantly between field-strength groups. In rat phantom data, differences between field strengths were less pronounced. At 7 T, substantial dataset-specific differences were observed in anatomical and functional quality metrics, also among datasets using the similar coil configuration. Multicenter reference values for frequently represented 7 T configurations were established, including SNR/tSNR of 48 {+/-} 3/48 {+/-} 1 dB for mouse surface-coil datasets and 49 {+/-} 3/44 {+/-} 4 dB for rat array-coil datasets. Longitudinal anatomical SNR showed generally low variability, whereas ghosting was more variable across time. Automated outlier detection additionally identified transient acquisition instabilities that were not always apparent by visual inspection. Together, these findings show that field strength and nominal coil configuration alone do not adequately characterize MRI performance and support standardized, longitudinal phantom-based QC combined with automated analysis as a practical approach for improving the reliability, transparency, and interoperability of multicenter preclinical MRI data.

Subject motion and gradient heating-induced frequency and phase offsets result in spectral misalignment during single-voxel MR spectroscopy (SVS) acquisitions. Several methods have been presented to align the spectra, but their performance on a 7T system is unclear. This study aimed to evaluate the practical importance of four retrospective correction methods, namely creatine fitting, residual water, spectral registration, and cross-correlation. SVS data were collected from 127 participants (88/39 female/male, 39{+/-}15 years) using a semi-localization by adiabatic selective refocusing sequence at 7T. Changes in the spectral linewidth, signal-to-noise ratio (SNR), similarity (mean similarity matrix (SImean)), and metabolite quantification (concentration, relative Cramer-Rao lower bounds (rCRLB)) after correction were evaluated. A Wilcoxon signed rank test was used to evaluate these changes. The p-values were adjusted using the false discovery rate (q<0.05) for multiple comparisons. For significant changes, the effect size was further calculated using the Rosenthal formula. The results demonstrated that all correction methods resulted in a significant decrease in the spectral linewidth and significant increases in the spectral SNR and SImean. The mean changes by using the four methods were 0.76-0.99 Hz for the spectral linewidth change (effect sizes:0.73-0.86), 5.79-6.95 for the SNR change (effect size:0.60-0.72), and 0.03-0.03 for the SImean change (effect size:0.87-0.87). The correction also significantly increased the estimates of the metabolite concentrations of PCr, Gln, Glu, GPC, Ins, NAA, Scyllo, and Tau (effect size:0.25-0.57), and significantly decreased the rCRLBs (effect size:0.26-0.42) of Scyllo, Tau, and Gly. These results indicate that each method provides practical benefits in improving the spectral linewidth, SNR, and similarity by showing a strong to very strong effect of the improvements (effect size>0.60). The correction only showed a weak to moderate effect of the changes in the metabolite quantification (effect size<0.60), suggesting that whether the correction results in more reliable metabolite quantification requires further validation.

Purpose: Proton density (PD) and macromolecular proton fraction (MPF) are fundamental MR parameters characterizing water and macromolecular content in tissues. Both parameters can be mapped using a time-efficient MRI protocol, though PD maps require additional corrections for the reception profile (RP) non-uniformity and global scaling. Existing PD mapping methods suffer from uncertainties related to PD normalization to reference tissue. This study aimed to develop a data-driven RP correction method for PD mapping based on MPF as an internal scale. Theory and Methods: One-parameter model linking PD and MPF was derived under assumption of proportionality between total macromolecular volume and proton concentration. This model with a uniform parameter value across brain tissues generates a surrogate PD map, which is used to estimate RP from an uncorrected PD map. The proposed algorithm was tested against a common PD mapping method utilizing bias-field correction and normalization to cerebrospinal fluid. Brain MRI data were acquired at 3T from seven healthy individuals. The protocol included 3D multi-echo T1-, PD-, and magnetization-transfer-weighted gradient-echo sequences. Parametric maps were reconstructed using the single-point MPF mapping algorithm with data-driven B1 correction. Results: MPF-based and bias-field corrections yielded similar PD measurements in white and gray matter with no significant bias and remarkably small variability (coefficients of variation <0.6%). Conclusions: MPF-based RP correction enables computationally simple and robust quantitative PD mapping at no additional cost if MPF maps are available. The main advantage of the proposed method is automatic elimination of global scaling that obviates the need for PD normalization.

Purpose: To evaluate, within a synthetic pulmonary-artery benchmark, whether biomarker-aware super-resolution improves hemodynamic endpoint fidelity in 4D flow MRI. Methods: A 3D residual channel attention network (RCAN) was trained on 22 CFD simulations from 11 pulmonary-artery geometries with 2x spatial degradation. The original primary comparison evaluated full biomarker-aware training against a combined divergence-temporal-regularization control across 14 matched seeds on a fixed six-case test set. Additional matched analyses subsequently separated divergence and adjacent-frame temporal-difference regularization and evaluated biomarker supervision in a controlled 2x2 design with a common reconstruction schedule. Results: In the original primary comparison, biomarker-aware training reduced PI error relative to the matched combined-regularization control (12.76 +/- 2.34% vs. 17.12 +/- 1.23%; paired difference -4.35 percentage points [95% CI -5.58, -3.09]; p = 0.0015). In the matched factorial analysis, biomarker supervision reduced PI error without (-2.235 percentage points [95% CI -3.514, -0.849]) and with (-4.125 percentage points [95% CI -5.392, -2.732]) combined regularization. PI improvement was not mirrored by uniformly improved reconstruction metrics: PSNR changed by -0.058 (p = 0.50) and +0.144 dB (p = 0.28), while SSIM decreased by 0.009 (p = 0.004) and 0.007 (p = 0.009), respectively. Divergence alone changed PI by +0.034 percentage points [95% CI -0.705, 0.809] while reducing DivRMS by 1.869 percentage points; temporal-difference regularization increased PI error by 1.562 percentage points [95% CI 0.594, 2.439]. Conclusion: Within this synthetic, fixed-case RCAN benchmark, biomarker supervision improved PI under both regularization conditions. PI improvement was not mirrored by uniformly improved reconstruction metrics. Temporal-difference - not divergence - regularization was the only isolated component associated with reproducible PI deterioration. These findings are configuration-specific and require geometry-independent and in vivo validation.

Purpose: Magnetic Resonance Acoustic Radiation Force Imaging (MR-ARFI) encodes focused ultrasound (FUS)-induced tissue displacements to MR phase shifts. Because these phase shifts are small, optimizing displacement sensitivity and minimizing motion artifacts requires characterization of the temporal dynamics of acoustic radiation force (ARF). Beyond sequence design, such characterization enables MR-ARFI to serve as a probe of brain viscoelastic behavior and can inform the design of therapeutic ultrasound pulsing schemes by quantifying how rapidly tissue deforms and relaxes. Methods: We developed a dynamic model assuming a gamma-variate impulse response for ARF(t), with time constant ({tau}) to be estimated. A single-shot spin echo ARFI acquisition with a bipolar motion encoding gradient (MEG) was implemented, incorporating retrospective corrections for head motion and physiological noise. A 500 kHz FUS pulse was applied in on/off blocks while systematically varying the timing between the MEG and sonication. The resulting phase time series were fit to the model to estimate time constant in five human participants. To assess sensitivity to tissue stiffness, {tau} was additionally measured in two customized gelatin phantoms of differing stiffness, characterized independently by MR elastography (MRE). Results: The average measured time constant in human participants was {tau} 3.7 {+/-} 0.27 ms. Comparing the two phantoms, the stiffer phantom exhibited a shorter {tau} (3.0-3.5 ms) than the softer phantom (4.0-4.5 ms), which tracked the MRE-derived Kelvin-Voigt viscoelastic time constants. In addition, the peak ARFI displacement was lower in the stiffer phantom (0.17 vs 0.74 m for the less stiff phantom), and a broader displacement focal spot size (FWHM 22 mm vs 15 mm). Conclusion: This time series approach with motion and physiological noise corrections validates an MR-based technique capable of characterizing the brain's temporal response to sonication-induced mechanical forces, providing a foundation for optimizing the design of ultrasound pulsing schemes for ARFI, as well as for future studies of brain biomechanics.

Purpose: To quantify, in a computational fluid dynamics (CFD) benchmark, how spatial degradation, temporal acquisition and noise change the error of plane-based pulmonary-artery flow endpoints. Methods: Twenty-two historical, non-converged CFD simulations (11 geometry families) were voxelised at 2 mm with 80 frames per cycle. Simulations with numerical divergence were excluded, leaving 13 (7 families). Fields were degraded to 4 and 8 mm and with noise, plane waveforms were sampled at 20 phases with one-interval averaging, and errors were computed against the noise-free reference waveform on 66 geometric planes. Family-mean paired differences are reported with t-based 95% confidence intervals under an analysis plan fixed in advance. Results: Spatial degradation from 2 to 4 mm increased full-waveform nRMSE by 0.0032 [0.0021, 0.0043], in 7 of 7 families; cycle-mean flow bias accounted for 23% of the squared error, against 60% at 8 mm. With SNR scaled to voxel volume, the increase was not resolved (0.0002 [-0.0012, 0.0015]). Twenty-phase acquisition increased PI error (0.0223 [0.0214, 0.0233]) and diastolic-trough timing error (0.63 [0.52, 0.74] frames) more than spatial degradation, and its point estimates for these endpoints were larger with every recovery method; the nRMSE ordering and peak timing depended on the recovery method. The implementation of the downsampling operator mattered: an alternative, corner-aligned implementation of the same 4 mm degradation gave a larger spatial effect (nRMSE difference 0.0010 [0.0002, 0.0017]). Conclusion: In this benchmark, a 2 to 4 mm degradation produced a small, consistent, mostly shape-related nRMSE increase that depended on the noise model, whereas 20-phase acquisition changed PI error and trough timing more. The results describe benchmark behaviour, not clinical acquisition requirements.

Purpose: to introduce a novel magnetic resonance imaging (MRI) technique, MUSLIS MRI (muscle and lipid suppression), that generates two co-registered 3D sets of images in a single acquisition: one set with lipid signal suppression (MUSLIS1) and a second set with muscle signal suppression (MUSLIS2). Materials and Methods: Numerical simulations based on Bloch equations were performed to optimize the sequence parameters; the MP2RAGE sequence was employed for the implementation of MUSLIS imaging. MRI experiments were performed in vivo at 3 T. Quantitative measurements of signal contrast were carried out on MUSLIS1 and MUSLIS2. Various combination images and T1 parametric maps of the thighs were generated from the two co-registered 3D sets. Results and Discussion: Optimal sequence parameters to achieve lipid suppression in one dataset (MUSLIS1) and muscle suppression in the other dataset (MUSLIS2) were determined by numerical simulations. MUSLIS imaging was validated in vivo at 3 T: MUSLIS1 images and MUSLIS2 images displayed excellent lipid and muscle suppression, respectively. The values of the signal contrast between subcutaneous lipids and muscle, as well as between bone marrow lipids and muscle reached values close to 0.8. Thus, MUSLIS would allow for tissue characterization of lesion/masses and direct visualization of the adipose tissue. Combination images, generated from these two intrinsically co-registered sets of images (MUSLIS1 and MUSLIS2), displayed different contrasts and improved signal characteristics, such as mitigation of signal inhomogeneities. The development of MUSLIS imaging for musculoskeletal imaging followed the concept introduced with the FLAWS, an MRI sequence designed for brain MR imaging to generate two 3D sets of images, one with white matter suppression and the other with cerebrospinal fluid suppression. Similarly to FLAWS, which provides T1 maps of the brain, MUSLIS imaging yields T1 parametric maps of the musculoskeletal system. Furthermore, the unique characteristics of the combination images make MUSLIS imaging a potentially attractive technique at field strengths higher than 3 T. Conclusion: The proof-of-concept of a novel approach to musculoskeletal imaging, MUSLIS MRI, is proposed here and demonstrated by numerical simulations and MRI experiments performed in vivo at 3 T. Keywords: MRI, muscle, lipids, musculoskeletal imaging, high field MRI, fat suppression, FLAWS MRI, MP2RAGE, MUSLIS, inversion recovery, adipose tissue, bone marrow, subcutaneous fat, T1 mapping, 3 T, 7 T, neuromuscular disorders, metabolic diseases, muscular dystrophy.

Purpose: To develop and evaluate a recurrent neural network that operates directly on undersampled golden-angle radial k-space, enabling low-latency cardiac cine reconstruction for interventional cardiac MRI (iCMR). Methods: KOSTAL (K-space Optimized Spatial Temporal Architecture for Low latency) combines a per-coil k-space convolutional gated recurrent unit (convGRU), which learns a mask-conditioned view-sharing rule over the accumulated k-space, with an image-space recurrent module. It was trained and tested on the open-access multi-coil OCMR dataset using 10 simulated golden-angle spokes per frame against a 1000-spoke reference. Ablations varied recurrent module architecture and placement; temporal behavior was assessed using a imaging-plane transition test. Generalization was evaluated without retraining in 23 pediatric patients prospectively scanned during free-breathing using a radial golden-angle bSSFP sequence with binned XD-GRASP and clinical breath-held cine serving as references. Results: At 10 spokes per frame, KOSTAL reached a median SSIM of 0.907 (gridding: 0.622) in under 34 ms of processing time per frame. Image quality (IQ) matched gridded results formed from 7-11 times more spokes. Recurrence in k-space rather than image space accounted for most of this gain (0.891 versus [&le;]0.650). IQ recovered from an abrupt plane transition in 75 ms, versus 303 ms for 200-spoke gridding. KOSTAL had the best agreement with the binned reference and the highest long-axis CNR. KOSTAL accurately tracked left-ventricular blood-pool volume (end-diastolic volume bias -0.7 mL) with a shorter lag (71 ms) than iterative reconstruction (239-287 ms). Conclusion: Learned, mask-conditioned recurrence in k-space yields high image quality at a latency compatible with interventional guidance.

PurposeWith growing interest in pediatric MRS to understand the neurochemistry of brain development, it is essential to establish pediatric-specific metabolite T1 and T2 relaxation times for reliable quantification of metabolites. Furthermore, characterization of age-related changes in relaxation times may provide insight into developmental differences in tissue microstructure and cellular composition. Therefore, the goals of this study were two-fold: (1) to estimate metabolite T1 and T2 relaxation times in a pediatric cohort; and (2) to assess linear relationships between age and relaxation times. MethodsTI- and TE-series data were acquired from the bilateral thalamus in 20 and 25 participants aged between 0 and 7 years, respectively. Data were processed and modelled in Osprey 2.5.0. Metabolite T1 and T2 relaxation times of total N-acetylaspartate at 2.0 ppm (tNAA2.0), total creatine at 3.0 ppm (tCr3.0), total creatine at 3.9 ppm (tCr3.9), total choline (tCho), myo-inositol (mI) and the combined glutamate and glutamine (Glx) were estimated by fitting the modelled signal amplitudes to TI-series inversion recovery signal and TE-series exponential decay equations, respectively. Linear relationships with age were tested for metabolite T1 and T2 relaxation times using Spearman correlations. ResultstNAA2.0 and Glx showed significant changes in T1 with age. tNAA2.0, tCr3.0, tCho and mI showed significant decrease in T2 with age. ConclusionSeveral metabolites showed significant changes in relaxation times with age in the thalamus region. Therefore, use of age-specific relaxation correction for pediatric studies using the linear equations provided is encouraged for reliable quantification of metabolites.

Patients with implantable cardioverter-defibrillators (ICDs) often are excluded from receiving CMR scans at 3.0T due to ICD-induced image artifacts. This work aims to transfer sequence protocols established for CMR imaging in the presence of ICDs at 1.5T to a first human application at 3.0T. All measurements were performed on a clinical 3T MRI scanner. Phantom measurements were used to establish and optimize wideband PSIR LGE imaging at 3.0T. Volunteer measurements with externally placed ICDs were used for optimization of CINE imaging and wideband PSIR imaging. Thirdly, patients with implanted ICD were examined with the optimized protocols. Human measurements showed that SA stacks obtained with spoiled GRE after contrast agent showed a markedly lower artifact burden than the conventional approach. For scar imaging, a wideband PSIR with a bandwidth of 3kHz and a frequency offset of +750Hz allowed a sufficient suppression of relevant artifacts in the majority of scans. No strong correlation between the image quality and the distance between the cardiac silhouette and the device was detected. The presented work shows the feasibility of clinical cardiac MRI at 3.0T including GRE-CINE and wideband PSIR-LGE for patients with ICDs.

Myelin supports rapid and efficient neural communication and is implicated in brain development, learning, aging, and neurological disease. Although numerous magnetic resonance imaging (MRI) approaches have been proposed to characterize myelin in vivo, they have rarely been evaluated together within the same participants and across repeated acquisitions. To address these gaps, we introduce the Myelin Imaging Reliability and ReprOducibility Resource (MIRROR), an openly available multimodal test--retest MRI dataset of diffusion MRI (dMRI), inhomogeneous magnetization transfer imaging (ihMT), quantitative susceptibility mapping (QSM), MP2RAGE, T1w/T2w ratio mapping, and transverse relaxometry in 22 healthy adults. We evaluated gray--white matter differentiation, intermetric correspondence, test--retest reliability, and participant specificity across 27 myelin-sensitive metrics. Gray--white matter differentiation was strong for dMRI derivatives, MP2RAGE-derived R1, T1w/T2w, and ihMT, while being substantially weaker for QSM measures. Intermetric spatial correlations were generally strongest within acquisition or modeling families. Test--retest reliability was generally high, particularly for dMRI-derived metrics, but was more variable among QSM-, ihMT-, and R2*-derived measures. Regional spatial profiles were also highly participant-specific, with discriminability approaching its upper bound for most metrics. All imaging data, derived maps, and processing code are openly released. MIRROR provides a benchmark for characterizing the measurement properties and practical tradeoffs of scalable myelin-sensitive MRI methods, supporting principled acquisition selection and accelerating methodological development and validation.

In magnetic particle imaging (MPI), signal is generated from the non-linear magnetic response of superparamagnetic iron oxide (SPIO) nanoparticles. The optimization of SPIOs for MPI is an active area of investigation. Tracer performance can be evaluated by magnetic particle relaxometry (MPR) and image-based metrics. This study characterized and compared the performance of ten commercially available tracers, spanning a broad range of magnetic core sizes, hydrodynamic diameters and surface coatings, using the peak signal intensity measured from MPR, and the total signal and maximum signal intensity measured from 2D MPI images. Synomag-D tracers exhibited the highest MPR peak signal intensities and highest maximum signal per {micro}g iron. ProMag, a micron-sized iron oxide particle (MPIO), yielded a low MPR peak signal intensity but had the highest total MPI signal per {micro}g of iron. A strong correlation was observed between MPR peak signal intensity and MPI maximum signal intensity (R2=0.94). In contrast, there was a weak correlation between MPR peak signal intensity and MPI total signal intensity (R2=0.38), though this correlation improved when MPIOs were excluded from the analysis (R2=0.91). While MPR generally predicts the tracer performance, it does not completely replicate the complex imaging conditions. As such, comprehensive tracer evaluation requires combining both MPR and image-based metrics.

Magnetic resonance imaging (MRI) exhibits substantial scanner and protocol dependent appearance variations, making harmonization challenging without distorting patient specific anatomy, particularly in abdominal imaging. We propose AC-Diff, an anatomy--contrast disentangled diffusion framework that formulates MRI harmonization as a factor-specific generative intervention. Using aligned multi-contrast supervision and cross-patient factor swapping, AC-Diff learns to separate anatomical structure from acquisition-dependent contrast. Unlike conventional image or latent diffusion, AC-Diff restricts stochastic generation to the contrast subspace while the source anatomy bypasses diffusion and is directly reused during reconstruction. Experiments on an in-house multi-vendor cohort and the external Duke Liver MRI dataset demonstrate improved target-domain alignment with strong structural preservation. In downstream liver segmentation, AC-Diff improves Dice from 0.942 for the original inputs to 0.954 for the harmonized images. These results support contrast-specific latent generation as a promising approach to anatomy-preserving MRI harmonization.

Automatic segmentation of medical images could accelerate clinical and research workflows, yet existing tools are rarely compared independently on shared benchmarks, making it difficult to assess genuine progress. This is particularly consequential for limb muscle segmentation in magnetic resonance imaging (MRI), which underpins volumetric analysis, radiomics, and fat fraction quantification used as biomarkers in neuromuscular disease -- applications where segmentation errors can directly affect interpretation. Despite the availability of multiple open-source segmentation tools spanning diverse architectures, no independent head-to-head comparison has appeared in the literature. We address this gap by evaluating eight tools, six of which run independently, and the other two of which are meant to be auxiliary algorithms once muscles have been somehow initially labeled. We evaluated these tools on three-dimensional MRI volumes across three base datasets and then later against derived datasets, assessing both quantitative metrics and qualitative usability. Our results reveal substantial performance variation across methods and datasets, with many methods showing reduced accuracy on pathological cases. Domain-specific models trained on large datasets consistently outperformed foundation models and newer general-purpose architectures. Our results reveal practical trade-offs between accuracy, generalizability, and ease of use that are critical for clinical adoption, while raising concerns for deployment on underrepresented populations and pathologies.

Quantitative magnetic resonance imaging (qMRI) supports both basic neuroscience and clinical research, providing the foundation for longitudinal studies of brain development, aging, and neurological disease. However, its validity rests on a little-studied assumption: that physical relocation of imaging hardware leaves underlying neuroimaging biomarkers intact. When research facilities expand or modernize, an entire scanner undergoes radical disruption--field ramp-down, physical transport, shielding reconstruction, and re-shimming. Whether qMRI metrics remain stable following such invasive procedures without introducing systematic biases is largely unquantified. Here, we systematically evaluated the impact of relocating a 3T clinical scanner over a distance of 30 km on a broad, multimodal suite of neuroimaging markers in healthy volunteers. To disentangle variance components, we benchmarked within- and cross-site stability using percent variation (PV; measurement precision), intraclass correlation coefficients (ICC; absolute reliability), and Spearman rank correlations (R2; spatial and order agreement). We show that scanner relocation introduces no greater variance than baseline within-site assessments. High-resolution brain structural morphometry from different sequences (multi-echo MPRAGE, compressed sensing MPR2RAGE) and quantitative relaxometry (compressed sensing MP2RAGE T1 mapping) demonstrated excellent stability (PV<1%, ICC=0.89-0.99, R2>0.93). Diffusion tensor (mean diffusivity: MD, fractional anisotropy: FA) and mean signal kurtosis (MSK) metrics from multi-shell diffusion MRI maintained low variability across white and deep gray matter (global cross-site ICC>0.90; global PV: MD<1.7%, MSK<1.2%, FA<3.8%; global R2: MSK>0.93; MD>0.90; FA>0.77), confirming relocation robustness for diffusion-derived microstructural markers. Resting-state functional networks displayed high cross-site reproducibility (ICC>0.90), while task-based fMRI activation peaks retained acceptable spatial and magnitude agreement (ICC=0.30-0.60) despite intrinsic physiological variability. In the absence of software or hardware upgrades, physical scanner relocation leaves qMRI biomarkers remarkably intact--providing the empirical foundation required to preserve data integrity across major facility transitions and multi-site longitudinal cohorts. Key points that summarize the paperO_LIHardware transitions do not compromise qMRI data integrity: Physical relocation of a 3T scanner (including field ramp-down, 30 km transport, and re-shimming) introduces no greater variance in quantitative MRI metrics than baseline within-site test-retest assessments. C_LIO_LIMultimodal biomarker robustness: Brain morphometry, T1 relaxometry, diffusion tensor and kurtosis, resting-state functional connectivity display comparable within- and cross-site stability reflected by different metrics (percent variation, intraclass correlation coefficient and Spearman correlation) C_LIO_LIEmpirical backing for longitudinal research: In the absence of hardware or software upgrades, multimodal qMRI biomarkers remain remarkably stable, reassuring researchers that long-term cohort studies and multi-site harmonization can survive major facility relocations. C_LI Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=175 SRC="FIGDIR/small/26363342v1_ufig1.gif" ALT="Figure 1"> View larger version (57K): org.highwire.dtl.DTLVardef@773d03org.highwire.dtl.DTLVardef@46a09borg.highwire.dtl.DTLVardef@da5bc0org.highwire.dtl.DTLVardef@1bdbf9d_HPS_FORMAT_FIGEXP M_FIG This study evaluates the impact of physically relocating a 3T MRI scanner on neuroimaging biomarkers. By comparing before and after relocation sessions on healthy volunteers, we conclude that the procedure does not introduce significant variability. Parameters for morphometry, T1 relaxometry, diffusion, and resting-state functional connectivity show high stability. Task-based fMRI showed greater within-site variability, but no increase in variability following relocation, consistent with the inherently higher variability of task-evoked responses.. The results confirm the integrity of qMRI data after infrastructure relocation, ensuring reliability for longitudinal studies. C_FIG

Early human cortical development involves rapid changes in cytoarchitecture that remain difficult to interrogate non-invasively. Biophysical diffusion MRI models provide a potential window onto this microstructure, including water exchange across cell membranes, but their assumptions are difficult to validate directly in developing human tissue. Human cortical organoids provide a tractable biological model in which diffusion measurements can be related to tissue architecture within the same specimen. We scanned six fixed cortical organoids at 9.4~T, three grown under a directed and three under an undirected differentiation protocol, at three diffusion times and b-values up to 7000~s/mm$^2$. Center-to-periphery profiles of DKI-derived diffusion metrics were reproducible across the directed organoids but more variable across the undirected organoids, while nuclear staining revealed greater variability in the internal architecture of the undirected organoids. Unexpectedly, mean diffusivity increased with diffusion time in more than 80% of tissue voxels in both batches, opposite to previous organoid measurements. Across 12,348 simulated two-compartment substrates, this increase was reproduced only with permeable membranes and an extracellular diffusivity substantially exceeding the apparent intracellular diffusivity, a regime that challenges assumptions commonly used in current gray-matter exchange models. These results support cortical organoids as biologically realistic platforms for probing early cortical microstructure and testing diffusion models under controlled conditions, while remaining amenable to post-MRI histological validation.